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ABSTRACT 

We investigate the observability of cold accretion streams at redshift 3 via Lyman- 
alpha (Lya) emission and the feasibility of cold accretion as the main driver of Lya 
blobs (LABs). We run cosmological zoom simulations focusing on 3 halos spanning al- 
most two orders of magnitude in mass, roughly from 10 11 to 10 13 solar masses. We use 
a version of the Ramses code that includes radiative transfer of ultraviolet (UV) pho- 
tons, and we employ a refinement strategy that allows us to resolve accretion streams 
in their natural environment to an unprecedented level. For the first time in a simula- 
tion, we self-consistently model self-shielding in the cold streams from the cosmological 
UV background, which enables us to predict their temperatures, ionization states and 
Lya luminosities with improved accuracy. We find the efficiency of gravitational heat- 
ing in cold streams in a ~ 10 11 solar mass halo to be around 10-20% throughout most 
of the halo but reaching much higher values close to the center. As a result most of 
the Lya luminosity comes from gas which is concentrated at the central 20% of the 
halo radius, leading to Lya emission which is not extended. In more massive halos, of 
> 10 12 solar masses, cold accretion is complex and disrupted, and gravitational heat- 
ing does not happen as a steady process. Ignoring the factors of Lya scattering, local 
UV enhancement, and SNe feedback, the cold 'messy' accretion alone in these massive 
halos can produce LABs that largely agree with observations in terms of morphology, 
extent, and luminosity. Our simulations slightly and systematically over-predict LAB 
abundances, perhaps hinting that the interplay of these ignored factors may have a 
negative net effect on extent and luminosity. We predict that a factor of a few in- 
crease in sensitivity from current observational limits should unambiguously reveal 
continuum- free accretion streams around massive galaxies at z = 3. 

Key words: cosmology: theory, diffuse radiation, large-scale structure of Universe, 
methods: numerical, radiative transfer 



1 INTRODUCTION 

The last decade has seen a shift in the way galaxies are 



thought to have assembled. In the classic theory ( Rees & Os- 
triker||1977||Silk||1977| [White fe Rees||1978| ), galaxies collect 
their baryons via so-called hot mode accretion where diffuse 
gas symmetrically falls into dark matter (DM) halos and is 
shock-heated as it hits the gas residing in them. Depending 
on the mass of the halo, the gas may or may not eventually 
settle into the galaxy. However, it has become increasingly 
apparent through theoretical work and simulations that at 
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high redshift (z > 2), galaxies get their baryons primarily 
via accretion of relatively dense, cold (10 4 K) and pristine 
gas which penetrates in the form of streams through the dif- 
fuse shock- heated medium (Far dal et al.pOOl) |Birnboim & ; 
|Dekel|2003||Keres et al.|2005||Dekel fe Birnboim|2006| |Birn- 
boim et al. 20071 lOcvirk et al. 120081 IDekel et al. 2009; Brooks 



et al 120091 |van de Voort et aI.|20lib||Faucher-Giguere et al 
2011||van de Voort &; Schaye|2011| )7Simulations consistently 
show these streams to exist and peak in activity around red- 
shift 3, though it appears that their widths are still dictated 
mostly by resolution. 

The problem is that cold accretion streams have never 
been directly observed, though we are starting to see some 
hints, both in emission ( Rauch et al.|[20lT ) and absorption 
( |Ribaudo et al.||2011| ). 
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Is this lack of observational evidence consistent with the 
existence of cold accretion streams? Do we not observe them 
because they're not easily observable or simply because they 
don't exist? 

Faucher-Giguere & Keres (2011) showed that the 



streams are hard to detect directly via absorption due to 
their small covering factor and surrounding galactic winds 
that overwhelm their signature. Kimm et al. (2011) came 



to the same conclusion, adding that the low metallicity in 
streams 10 -3 solar) further inhibits their detection via 



metal line absorption. Even so, Fumagalli et al. (2011) and 
van de Voort et aL| ( 2011a| ) have argued that a large fraction 



of observed metal-poor Lyman-limit systems (LLSs) make 
up for indirect detections of cold streams. Furthermore, we 
may possibly have been directly observing the tips of these 
streams during the last decade in the form of Lyman-alpha 
blobs (LABs). 

LABs are extremely bright (> 10 43 ergs -1 ) and ex- 



tended (> 30 kpc in diameter) Lya nebulae (e.g. Francis 



eT"aL|[T996l |Keel et al.||1999| |Steidel et al.||2000l |Matsuda 
et al.|20041|Palunas et al.||2004| |Nilsson et al.||20061 |Smith| 
fe Jarvis||2007| |Prescott et al ||2009| |Yang et al.||2010| |E"rb| 
et al.||2011 ). They have a slight tendency to be filamentary 



in structure ( Matsuda et al.|20lT hereafter Mil), and often 
have short limbs protruding from the main body. They of- 
ten coincide with galactic sources that give hints about their 
physical origin but the mechanism by which the emission be- 
comes so strong and extended is a matter of debate. A sub- 
set of LABs however have no apparent coinciding galactic 
sources (e.g. |Steidel et aT1|2000| |Weijmans et al.||2010| |Erb| 
|et al.|[20lT ). Up until now about two hundred LABs have 
been discovered, including about fifteen giant ones (> 100 
kpc). Smaller extended Lya emitters exist in large quanti- 
ties over a continuous range of sizes down to point sources. 
LABs appear to be specific to the high-redshift Universe 
(Keel et al. 2009) and most of them have been detected at 
2 < z < 3. 

The physical nature of LABs is still a matter of debate, 
but by most accounts they are powered by a combination of 
some or all of the following processes: (a) Cold stream ac- 
cretion is a natural explanation, where the fuel source is the 
dissipation of gravitational potential, also termed gravita- 
tional heating (e.g. |Steidel et al.|2000| |Haiman et al.|2000| 



Fardal et aT1|200T] [Dijkstra et al.||2006| |Dijkstra Loeb| 
2009). (b) Photo-fluorescence by near- lying sources, such as 



active galactic nuclei (AGN) or starbursts (e.g. Haiman & 



Rees|2001|[Cantalupo et al.| 2005; K ollmeier et al.|2010| ), Jcj 

Lya scattering, also fuelled by neighbouring star-forming 
regions (e.g. Laursen h Sommer-Larsen 2007 Zheng et al. 
2011). (d) Cooling radiation in galactic outflows, fuelled by 
AGN or supernovae (e.g. Taniguchi &; Shioya 2000; Ohyama 



et al.||2003||Mori et al.|2004| ). 



Furlanetto et al. (2005) used cosmological simulations 



to look at the contributions of each of these processes, and 
found that star-forming regions can in principle power all 
but the largest LABs via photo-fluorescence and Lya scat- 
tering, but that cold accretion alone cannot, except under 
very optimistic assumptions. They however pointed out that 
the Lya emissivity of their simulated gas is highly uncertain 
due to the lack of modelling of self-shielding from UV ra- 
diation: The self-shielding state of the gas affects both the 



temperature and ionization state, which sensitively dictates 
the Lya emissivity. They also pointed out that the efficiency 
of star-formation in powering LABs is very dependent on 
the presence of dust. As pointed out by Cen (2011 ), massive 



galaxies tend to have large dust content which makes them 
very efficiently transform their UV (and Lya) output into 
infrared radiation. Thus it appears problematic to associate 
the largest and most luminous LABs to star-formation in 
the most massive halos in the Universe. 



1.1 Recent work on gravitationally driven Lya 
emission 

Notably, two recent simulation papers have studied gravi- 
tational heating as the driver of LABs, but have reached 
conflicting conclusions: 

Goerdt et al. (|2010[ hereafter G10) analyze two suites 



of cosmological adaptive mesh refinement (AMR) simula- 
tions. They assume self-shielding in post-processing from 
the UV background in accretion streams. Mock observa- 
tions of halos of rsj 4 10 n MQ at redshift 2.3 look similar 
to real LABs in morphology and surface brightness profile, 
though the association of LABs to halos of such low mass 
implies an unrealistically high LAB abundance. A Lya lumi- 
nosity function derived from their results is not far off from a 
function derived from observations, though they over-predict 
number densities somewhat, which implies the cooling emis- 
sion in their simulations is too efficient. As pointed out by 
Faucher-Giguere et al. ( 2010[ hereafter FG10) this overesti- 
mate appears to be due to artificial photo-heating of stream 
gas, which is not on-the-fly self-shielded from the UV back- 
ground. 

FG10 analyze cosmological smoothed particle hydrody- 
namics simulations to test different approaches and approx- 
imations. Based on radiative transfer (RT) post-processing 
results, they apply on-the-fly self-shielding by excluding UV 
photoionization from all gas denser than 10 ~ 2 H atoms per 
cm 3 . Then they apply a Lya transfer code to their output to 
model the scattering of Lya photons towards the observer 
and obtain realistic mock observations. According to their 
results, which are in good agreement with Furlan etto et al.| 
(2005), cooling radiation can in principle power LABs, pro- 
vided one includes emission from gas dense enough to be 
star-forming to some extent. They note that this gas should 
be under the influence of feedback processes which introduce 
a large uncertainty to the cooling emission. 

Although G10 and FG10 are not in good agreement 
on their conclusions, they both agree with Furlan etto et al.| 
fl2005| ) on that proper modelling of self-shielding from UV 
radiation is crucial to the results. 



1.2 This work 

We have developed a radiation- hydrodynamics (RHD) ver- 
sion of the AMR code Ramses (Teyssier 2002), which puts us 
in a unique position to continue the work of the aforemen- 
tioned authors, to study the emissivity of accretion streams 
in their natural environment at high redshift in simulations 
that accurately and consistently model self-shielding from 
the UV background. We also extend previous work by sim- 
ulating halos of larger masses, which are more likely to host 
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LABs, and by using an original refinement strategy which 
allows us to describe cold streams with unprecedented reso- 
lution. The increased resolution also allows us to accurately 
track the state of the gas up to higher densities than in 
the previous works. The main motivations of our work are: 
(a) Investigate whether gravitational heating is capable and 
sufficient as a driver of observed LABs. (b) Predict the ob- 
servability of gravitationally powered Lya emission from ac- 
cretion streams at redshift 3. 

The paper is structured as follows: Section 2 describes 
the simulation code and the setup of our experiments. Sec- 
tion 3 describes the physical properties at redshift 3 of our 
simulated halos over a range of masses. Section 4 presents 
our prediction of the Lya emission from extended gas around 
galaxies and its observability. We compare with observations 
of LABs. In section 5 we discuss the efficiency of gravita- 
tional heating as a source of extended Lya emission and the 
contribution of cosmological UV fluorescence. We discuss 
other factors that may affect the extended Lya emission. 
Finally we conclude and discuss in section 6. 



2 SIMULATIONS 
2.1 Code details 

We run our simulations in RamsesRT, a version of the AMR 



code Ramses (Teyssier 2002) which we have modified to 



include on-the-fly radiation-hydrodynamics describing the 
propagation in space of UV photons and their interaction 
with gas via photoionization and heating of hydrogen and 
helium. 

The widely used Ramses code simulates the cosmological 
evolution and interaction of dark matter, stellar populations 
and baryonic gas, via gravity, hydrodynamics and radiative 
cooling. The gas evolution is computed using a second order 
Godunov scheme for the Euler equations, while trajectories 
of collisionless DM and stellar particles are computed using 
a Particle-Mesh solver. 

The RamsesRT implementation and tests will be fully 
described in Rosdahl et al. (2012, in preparation), and here, 
we only briefly present the aspects of RamsesRT which are 
most relevant to the present work. 

For the radiative transfer we use a moment-based 
method with the Ml closure relation, as described in[Aubert 
& Teyssier (2008), the essence of which is to turn rays of 
radiation into a fluid with a direction of flow that corre- 
sponds to an average of rays over all angles. In contrast to 
the usual ray-tracing codes currently on the market, this 
gives the advantage that the computational load of RT does 
not scale linearly - in fact hardly scales at all - with the 
number of radiative sources in the simulation. This is a par- 
ticular advantage here as we simulate a spatially continuous 
source of radiation, which is hard to do with a ray-tracing 
code. RamsesRT takes full advantage of the AMR structure of 
Ramses and photons are propagated through the same cells 
that define the baryonic gas. 

Our RT solver is explicit, which means the timestep 
length for the propagation of photons is limited by the speed 
of light. This typically makes the RT timestep three orders of 
magnitude shorter than the hydrodynamical timestep. Since 
we're forced to apply a global RHD timestep which is the 



minimum of the hydrodynamical step and the RT one, we're 
faced with the rather horrifying prospect that RamsesRT sim- 
ulations are slowed down by a factor of order one-thousand 
compared to non-RT simulations. To get around this, we 
invoke the Reduced Speed- of- Light Approximation (RSLA) 
proposed by |Gnedin &; Abel| ( 2001| (see also discussion in 
|Aubert &; Tey ssier 2008 ) : The speed of light is reduced by a 
factor f c , bringing the RHD timestep closer to the normal 
Ramses one and making RamsesRT runnable in reasonable 
time. In the simulations described here we use f c ~ 1/100. 
To be sure the choice of light speed is not affecting our re- 
sults, we have run analogs of our HI simulation (see Table 
[TJ with the light speed changed by a factor of five in ei- 
ther direction, i.e. f c = 1/20 and f c — 1/500. This has an 
insignificant effect on the results, and we conclude it is an 
acceptable approximation for our simulations. 

In order to self-consistently evolve the UV field we im- 
plement non-equilibrium gas cooling that keeps track of the 
abundances of all ion species of hydrogen and helium. These 
abundances are stored in the form of three ionization frac- 
tions, as passive scalars that are advected with the gas, 



xmi 

XUell 
XHelU 



nmi/nu, 

^Helll/^He, 



(i) 



where n is number density. The non-equilibrium cooling 
module evolves these ionization fractions along with pho- 
ton fluxes and temperature on a per-cell basis, with the 
timestep constraint that none of these quantities changes by 
more than 10% in a single timestep, using sub-cycles when 
needed to fill the RHD timestep. 

We have tested and verified RamsesRT with the bench- 
mark tests of the 'Cosmological radiative transfer compari- 
son project' ( Iliev et al.|2006| 2009| ), and the results will be 
presented in Rosdahl et al. (2012, in prep.). 



2.2 Simulation setup 

We run three cosmological zoom simulations, each targeting 
the evolution until redshift 3 of a single halo and its large- 
scale environment. The initial conditions are generated us- 
ing MPGRAFIC flPrunet et al.|2008[ ). We assume a ACDM 
Universe with Q A = 0.723, Q m = 0.277, Vt h = 0.0459, 
h = ifo/100 = 0.702 an d a 8 0.817, consi stent with seven- 



year WMAP results ( Komatsu et al.[2011 ). We assume hy- 
drogen and helium mass fractions X = 0.76 and Y = 0.24. 

The masses of these three halos span almost two orders 
of magnitude, the least massive halo (~ 3 10 11 Mq) roughly 
corresponding to halos studied in G10 and FG10, and the 
more massive halo simulations (up to mass ~ 1O 13 M ) 
based on the expectation that LABs are situated in over- 



dense regions of the Universe (Steidel et al.||2000 Prescott 



et al. 2008 ; Y ang et al.|2010 ). The parameters for the individ- 
ual simulations, named HI, H2 and H3 in order of increasing 
halo mass, are listed in Table ^ 

Each simulation has periodic boundaries and nested lev- 
els of refinement in a zoom-region around the targeted halo, 
both in DM and gas. 

On-the-fly refinement is enforced inside the zoom re- 
gions according to two criteria: The first is the traditional 
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Table 1. Simulation parameters 



Name 


Box 
siz<0 

[Mpc] 


Halo 
masq^j 

[M ] 


DM 
res [£] 

[M ] 


Gas 

res0 

[pc] 


m 


[cm - " 


HI 


28.5 


2.9 10 11 


1.4 10 6 


217 


1/100 


i io- 4 


H2 


28.5 


2.9 10 12 


1.1 10 7 


434 


1/300 


3 IO" 4 


H3 


51.2 


1.3 10 13 


6.4 10 7 


780 


1/300 


3 IO" 4 



a Co-moving 

DM+baryons at z 



3, (all the mass within the virial radius) 

Optimal resolution 
Optimal physical resolution (not co-moving) at z = 3 



c 
d 

e Reduced light-speed fraction, see Sec 
f Threshold for UV-emitting gas, see Sec. 



2.1 



2.2 



' quasi- Lagrangian' criterion, where a cell is refined if it con- 
tains more than 8 DM particles or the equivalent baryonic 
mas£] This causes concentrations of mass to be refined to 
the maximum, but will typically leave the resolution of cold 
flows many times less, which is a problem when one is most 
interested in these flows. The second refinement criterion, 
which is unique to this work, is applied on the hydrogen 
neutral fraction gradient. According to it, two adjacent cells 
at positions i and 2 + 1 are refined (up to the maximum level 
of refinement), if 



i _|_ *+l _|_ floor 

x m -t- x m -+- x m 



> Axhi, 



(2) 



where xhi = 1 — ^hii and floor on the neutral 

fraction under which the criterion becomes inactive. In our 
simulations we typically use Axhi = 0.75 and x^° or — 10 -3 
in order to resolve gas streams, though we tweak these val- 
ues a bit (even within the same simulation) to tread the fine 
line of neither under-resolving the streams nor over-resolving 
uninteresting regions. This enforces maximum refinement in 
the cold streams, so while the optimal resolution in our sim- 
ulations is slightly less than in recent works, our resolution 
in the streams is unprecedented in cosmological simulations. 

The cosmological UV background is incorporated 
into our simulations with an 'outside-in' method, where it is 
propagated from under-dense and transparent UV-emitting 
voids. As such, our UV background can be thought of as 
quasi-homogeneous, as opposed to the completely homoge- 
neous and optically thin implementation commonly used 



in cosmological codes that lack radiative transfer (e.g. Cen 
~ Katz et al.|1996||Rasera fe Teyssier|2006| ) . The reasons 



1992 



we chose this model are mainly twofold. First, it is only a 
single step further than previous work on the subject. This 
allows us to isolate the effect of self-shielding, and to inter- 
pret our results in a well established theoretical framework. 
Second, the inside-out method would require finely tuned 
star formation rates and UV escape fractions for simulated 
galaxies to produce a 'correct' UV background, and this is 
a subject onto itself (see e.g. |Wise &; Cen||2009] Aubert & 
Teyssier|2010| . Also, our simulations zoom in on a relatively 



1 A cell is refined if it contains a mass of baryons larger than 
8 flfr/flm rriDMi where fl^ and Q m are the cosmological mass 
fractions of baryons and matter, respectively, and m^M is the 
mass of the highest-resolution DM particles. 



small volume with no star formation outside, which would 
lead to a severe lack of external UV background radiation. 
We thus postpone such a model to a future paper and in- 
stead demonstrate in Sec. 15.41 that a local enhancement of 
UV radiation due to star formation would not significantly 
change our conclusions. 

In practice, we use a 'void' density threshold n^ v such 
that all gas cells lower in density are UV emitters, and we 
impose the redshift-dependent UV background model from 
|Faucher-Giguere et al. ( 2009 ) onto these cells under the valid 
assumption that voids are optically thin. The radiative field 
is then allowed to diffuse out towards denser regions. The 
idea is to have the void threshold low enough that it doesn't 
include the potentially shielded cold streams themselves, but 
high enough that radiation can quickly reach them (a sen- 
sitive issue due to the reduced speed of light). We use void 
thresholds of n^j V > 10 -4 cm -3 in our simulations. Our re- 
sults are not sensitive to the fine-tuning of this as long as 



10" 



- 3 < 



i u H v < IO" 2 cm" 



The spectral shape of the UV field is approximately 
taken into account by using three (Hi-, Hei- and Heii- 
ionizing) packages of photons which are propagated indepen- 
dently (see Appendix [A| . In this work we adopt the on-the- 
spot approximation (OTSA), where any UV photon emitted 
from a recombination is assumed to re-ionize a nearby atom 
(i.e. within the same grid cell) - in other words, the simu- 
lated gas does not emit UV photons due to recombinations 
and case B recombination rates are used in computing the 
gas cooling rate. 

For the sake of simplicity, our simulations do not include 
SN feedback or metals. To prevent artificial fragmentation 
dTruelove et al.||1997 ) our simulations employ a polytropic 
equation of state ( |Dubois &; Teyssier] [2008 ) as a subgrid 
recipe that keeps the mostly unresolved multi-phase inter- 
stellar medium (ISM) from collapsing and fragmenting. The 
recipe sets a density-dependent temperature floor in every 
gas cell: 



To 



nun 



7-1 



(3) 



where we've chosen the values To = 10 K, nn,o = 1 cm - , 
and 7 = 1.6. The value of nu,o also corresponds to the limit 
above which gas is star-forming. 

We identify halos in our simulation outputs with the 
AdaptaHOP algorithm from |Aubert et al.| ( [2004| and |Tweed 
et al. ( 2009| , where the virial radius of a halo, R V ir, is de- 
fined as the radius where the average density is 200 times 
the critical density of the Universe, and the halo center cor- 
responds to the DM density maximum. 



2.3 Numerical issues 

We've established through convergence tests that resolu- 
tion is adequate in our simulations and that the chosen pa- 
rameters of light speed and UV emission threshold (f c and 
n H V ) do not affect our results noticeably. Three other issues 
should be noted: 

The gravitational potential in our simulations is 
usually dominated by DM particles, but it is resolved to the 
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local cell resolution. With our strategy of optimally resolv- 
ing gas streams comes the danger that we may over-resolve 
the gravitational potential, with few and far-between DM 
particles causing discreteness effects in the potential, which 
may lead to artificial fragmentation and complexity in the 
streams. This seems particularly ominous since we find in 
our simulations that the streams indeed become fragmented 
and complex around massive halos. To make sure this is not 
caused by an over-resolved gravitational potential we have 
run analogues to our simulations with smoothed potentials, 
which still reveal fragmented streams. So while it is hard to 
tell whether or not this numerical effect is nonexistent in our 
simulations, we can conclude that it is not dominating our 
results and that the complex streams are physical in nature. 

Operator splitting is a widely-used method of decom- 
posing unwieldy differential equations into separate parts 



that can be solved independently and in sequence (e.g. Toro 



1999 



Press et al.||1992| . RamsesRT employs this method to 



split the radiative- hydrodynamics equations into (i) advec- 
tion of gas and photons between cells and (ii) chemical reac- 
tions within the cells (radiative cooling and photo-heating) . 
The advection part is first solved and then cooling, using 
the advection result as the initial state and subcycling when 
needed. Gas normally exists in a competition between ad- 
vective/gravitational heating and radiative cooling, where 
the temperature 'adjusts' to a value where these processes 
cancel each other out. However, when the cooling time is 
shorter than the advection time, operator splitting may ar- 
tificially give cooling the upper hand, leading to a slight 
underestimate in the temperature. Normally this is no big 
deal, but considering how sensitive Lya emissivity is to gas 
temperature (see Fig. [6]), this can result in a severe underes- 
timate of Lya emissivity in the gas. We've verified that this 
is indeed the case in our simulations. To prevent this from 
affecting our results, we restart the simulations at z — 3 
with the global timestep reduced by orders of magnitude, to 
make sure it is everywhere shorter than the local cooling- 
time, and run until we reach convergence in Lya luminosity 
(this takes a few-thousand fine-cell time-steps). 

Cell merging: With the bookkeeping on ionization 
states, and due to the fact that cell de-refinement takes place 
just before outputs are written in Ramses, special care must 
be taken on cell merging. Applying the traditional method 
of giving a merged cell a children-averaged ionization state 
can sometimes result in a combination of temperature and 
ionized state which causes it to outshine whole galaxies in 



Lya emissivity (see discussion in Sec. 4.1). To prevent this 
we enforce a photoionization equilibrium (PIE) ionization 
state to merged cells, assuming the children- averaged values 
of gas density, pressure, and UV flux. 



3 PHYSICAL PROPERTIES OF 3 HALOS 

In this section, we first review the qualitative properties of 
our three simulated halos, and define the different phases of 
the intra-halo gas. We then describe in detail the impact of 
self-shielding on the ionization and thermal states of cold 
streams, and discuss the validity of an approximate treat- 
ment of self-shielding introduced by FG10. 



3.1 Basic halo properties 

Gas density maps of the three targeted halos at redshift 3 are 
shown in Fig. [I] in close-ups of the halos and zoom-outs to 
show their environments. Also shown are zoom-out maps of 
temperature. The halos display a tendency with increasing 
mass towards more intense, complex and fragmented accre- 
tion, and larger and hotter domains of shock-heated inter- 
galactic medium (IGM). 

The least massive halo (HI, left) has narrow (down to 
~ 1 kpc in diameter) and unperturbed accretion streams 
and tidal tails stretching from the central galaxy, which can 
be seen in red at the center of the halo, edge-on and slightly 
inclined from the horizontal. It is about to undergo a ma- 
jor merger with another halo five times less massive, situated 
just outside R V i r and coming in from the north. Two parallel 
accretion streams bridge the merging halos. Another accre- 
tion stream extends towards a factor 100 smaller merging 
halo, also at the edge of R V i r , but towards the line of sight 
(LOS), seen as a moon-shaped clump to the south-west. To 
the south and south-east are two relatively thick and diffuse 
accretion streams and another one even more diffuse to the 
west. Other structures in the map are orbiting satellites and 
tidal tails. 

The intermediate mass halo (H2, middle) is more a 
group of orbiting galaxies than a single galaxy. On the large 
scale there is a network of filaments mixed with galaxies of 
varying masses, with at least 6 large scale streams extend- 
ing towards the central halo. Movies show that the accretion 
here is notably more spiral than around the HI halo, with 
the streams starting to curve around the center of gravity 
already well outside R V i r . Inside the halo we see plenty of 
streams and tidal tails, but much more disrupted and messy 
than in the HI halo, as a result of stronger and more frequent 
interactions with other streams and galaxies. 

This tendency continues with the most massive halo 
(H3, right), where we find even more disrupted streams, to 
the point that many of them seem to be completely oblit- 
erated close to the halo center. The H3 halo has just un- 
dergone a major merger, which makes the accretion activity 
particularly violent at this point in time. 

To facilitate our analysis, we apply the following cat- 
egorization to divide the gas into phases, as shown in the 
temperature-density phase diagram in Fig. [5] (note that the 
categorization is specific to this paper and does not apply 
in general): 

The star-forming ISM is all gas denser than 1 cm -3 
and as discussed in Sec. |2.2| we apply a temperature floor 
in the form of a density-dependent polytrope to keep this 
gas from artificially fragmenting, which manifests itself in 
the constrained temperature-density relation in the shaded 
area of the diagram. Our simulations lack the ingredients to 
accurately model Lya emission from the ISM (multiphase 
resolution, dust, Lya scattering) and our reaction to that is 
to simply ignore the Lya emission from there in our analysis. 
The shaded color of the ISM region in Fig. [2] should remind 
the reader of this and that this work is about modelling the 
Lya emission coming from galactic environments and not 
the galaxies themselves. The ISM gas density threshold is 
resolution dependent and reflects the density at which fur- 
ther collapse of gas - i.e. the Jeans length - is no longer 



© 0000 RAS, MNRAS 000, 000-000 



6 Rosdahl & Blaizot 




HI H2 H3 



Figure 1. Redshift 3 maps of the three targeted halos in simulations HI, H2 and H3 from left to right (increasing halo mass). Grey 
circles indicate virial radii of the halos; 46, 98 and 158 kpc for the HI, H2 and H3 halos respectively. Top row: Number density maxima 
along the line of sight, with contours marking 0.02 and 0.3 cm -3 as indicated in the color bars, corresponding to our definition of streams. 
Middle row: The same but zoomed out to show the large-scale environment. Bottom row: Mass weighted temperature maps, on the 
same scale as the middle row. 

resolved. At our chosen density threshold, assuming mini- 
mum temperatures of 10 4 K, the Jeans length is resolved by 
approximately 10, 5, and 2.5 cell widths in the HI, H2 and 
H3 simulations respectively. It should be noted that our den- 
sity threshold is almost an order of magnitude above what 
has typically been used in recent similar works (e.g. G10, 
FG10). 

The CGM is gas with number densities between 0.3 



and 1 cm -3 . Ideally these densities form membrane inter- 
faces between the ISM and their more diffuse environment. 
The lower density limit corresponds to the inner contours in 
the density maps of Fig. [I] (top row) , and from those maps it 
can be seen that the CGM gas is indeed mostly constrained 
to galaxies (in red). In the phase diagram we find that most 
of the CGM gas is cooled down to the temperature floor of 
~ 10 4 Kelvin where radiative cooling basically stops (metals 
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Table 2. Halo sizes and mass budgets (% of mass within R v i r ). 



log(n H [cm" 3 ]) 

Figure 2. Phase diagram of the HI halo, showing our density- 
dependent definitions of streams, CGM and ISM. The ISM is 
shaded to indicate that we always ignore ISM gas when adding 
up Lya emissivities. The color scale represents mass weighted 
probability per temperature-density bin. 



can cool gas further but we don't include those). Although 
CGM gas in our simulations is not directly affected by the 
polytropic equation of state, one may expect gas at densi- 
ties > 0.1 cm -3 to be multiphase and star-forming JSchaye 



2004). This cannot happen in our simulations because they 
don't describe cooling below ~ 10 4 K, and this temperature 
floor provides artificial pressure support for dense gas. This 
implies a potentially high error in our predictions for the 
Lya luminosities of halos, resulting from an overestimated 
CGM contribution. Thus, while we in general include the 
CGM gas in our analysis of Lya luminosities, we also con- 
sider at some points the effect of excluding it, to get a grip 
on how sensitive our results are to the density thresholds 
applied. In summary we find that GCM gas typically pro- 
vides a 40% of the Lya luminosities of our simulated halos, 
but that in terms of Lya extent it is less substantial. 

Streams are denned in this work as gas with densities 
between 0.02 and 0.3 cm -3 . These limits correspond to the 
density contours in Fig. [I] (top row) and from those we can 
see that these densities indeed correspond to thin filamen- 
tary structures. Much like the CGM most of the stream gas 
is found at the bottom of the temperature curve at ~ 10 4 K 
though we do see an increase in temperature in the more dif- 
fuse stream gas due to a combination of photo/gravitational 
heating and inefficient cooling (because of the low densi- 
ties). Gas at sub-stream densities turns out to be negligible 
in terms of Lya emissivity and thus not very important to 
our results so we crudely split what remains into two cate- 
gories. 

Hot diffuse gas has been shock heated above 10 5 
Kelvin. As seen in the temperature maps in Fig. [I] (bottom 
row) this gas exists in abundance within the virial radii of 
the halos, but there also seems to be weaker heating around 
the large-scale accretion streams (and actually not so weak 
in the large streams around the H3 halo). Shock heating 
gets decidedly stronger with increasing halo mass, with gas 
reaching - 2 10 7 K in H2 and - 6 10 7 K in H3. Also, in- 
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creasingly dense gas exists above 10 K in the more massive 
halos; CGM in H2 and ISM in H3. 

Cold diffuse gas is partly gas which is slowly condens- 
ing towards the streams and the CGM and partly cosmolog- 
ical gas that has not interacted with the halos at all and is 
being cooled down by the cosmological expansion. 

The sizes and mass budgets of our three targeted halos 
are listed in Table |2j Each of the halo masses consists of 
roughly 80% dark matter and 20% baryons. The stellar/gas 
ratio decreases with halo mass, going from roughly one-to- 
one in HI to about a one-to-three in H3. The gas mass is 
primarily in the ISM, going from 73% in HI to 58% in H3. 
The hot gas fraction clearly increases with halo mass, going 
from 8% in HI to 23% in H3 and correspondingly the cold 
fraction decreases, going from 8% to less than 1%. Inter- 
estingly the stream fraction peaks in the intermediate mass 
halo at 16%, with half and two-thirds of that in HI and H3 
respectively. The low fraction in HI can be explained by the 
smooth accretion that efficiently moves the gas straight into 
the ISM (hence a high ISM fraction), whereas H3 streams are 
disrupted to the point of obliteration when they approach 
the halo center (hence the low ISM fraction and large frac- 
tion of hot gas) . 



3.2 On-the-fly self-shielding 

The transfer of UV photons gives us the opportunity to 
study the extent of self-shielding in gas clumps and streams. 
We quantify the local UV field intensity in terms of the 
hydrogen photoionization rate T, which expresses the aver- 
age number of photoionization events per hydrogen atom 
per unit time (see Appendix [5]). In the UV model we 
use, the unattenuated pho toio nization rate at redshift 3 is 
r = 6.1 10 -13 s _1 (see Fig. A2), and shielded regions should 
have r s" 1 . 



Fig. [3] shows the UV attenuation in the three targeted 
halos at redshift 3. The top row contains no n- logarithmic 
maps of projected minima of the photoionization rate along 
the LOS. The light color on the edges of the maps corre- 
sponds to the unattenuated value. Towards the centers of 
the halos the UV field becomes increasingly attenuated due 
to photo-absorption of the gas and in the densest streams 
and clumps we see ~ 100% attenuation. The diffuse streams, 
with densities < 0.02 cm -3 , are not self-shielded. Gas at the 
centers of the HI and H2 halos is efficiently shielded but at 
the center of the H3 halo gas is thermally ionized and thus 
optically thin. 

The bottom row of Fig. [3] shows logarithmic phase di- 
agrams of the hydrogen photoionization rate Y versus den- 
sity for the same three halos. The most diffuse gas is UV 
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Figure 3. Self-shielding at redshift 3 in the three halos, from left to right, HI, H2 and H3. Top row: Maps of projected minima along 
the LOS of the hydrogen photoionization rate. The scale is non-logarithmic and in units of 10 -13 s _1 . Bottom row: Phase diagrams 
of hydrogen photoionization rate versus density. The color scale represents mass-weighted probability per T — nu bin over the plotted Y 
range. 



emitting and has corresponding horizontal lines in the dia- 
grams up to the -threshold. Above this threshold there 
is an immediate spread in the photoionization rate in all 
three halos, ranging from unattenuated UV to about half 
attenuated. Gas in the HI halo is mostly self-shielded at 
riH > 2 10 -2 cm -3 . In more massive halos, the advent of 
thermal ionization in dense gas makes the situation more 
complex, and gas at uh ^ 2 10 -2 cm -3 exists in two 
phases, either self-shielded as in HI or optically thin. The 
bifurcation in the diagram around Y ~ 10~ 17 - 10~ 15 s~\ 
riH ~ 0.1 — 1 cm -3 , which grows more conspicuous with 
increasing halo mass is an effect of the ionization fronts be- 
coming under-resolved at high densities, where the mean 
free path becomes comparable or shorter than the cell sizes. 
This feature does not affect our results. 

Fig.[4]shows maps of Hi column densities and phase dia- 
grams of neutral fraction xhi = nm/nn vs. density. We find 
the CGM and ISM regions correspond mostly to damped 
Lya absorbers (DLAs, Nn > 2 10 20 cm -2 ) and the streams 
to Lyman limit systems (LLSs, Nu 1.6 10 17 - 10 19 cm" 2 ) 
and even super Lyman limit systems (SLLSs, Ah = 10 19 - 
2 10 20 cm -2 ), according to the definitions found in e.g. Fu 



magalli et al. (2011). The column densities are likely over- 



hanced UV from star-formation. We see in the phase dia- 
grams an abrupt transition of cold gas from ionized to neu- 
tral states, at about 5 10 



2 cm 3 in all halos. This generic re- 



sult is consistent with early expectati ons from|Schaye|(|2 001 ) 
and recent numerical estima tes (e.g. |Kollmeier et al.| |2010; 



estimated where they are highest due to lack of locally en- 



Faucher-Giguere et al ||2010| |Aubert fe Teyssier||2010| . The 

dense (tih ~ 0.1 cm - ) and ionized (xni < 10 _ ) cells which 
become increasingly abundant with halo mass correspond to 
hot shock-heated gas, which is thermally ionized and opti- 
cally thin. 



3.3 A self-shielding approximation 

FG10 applied a self-shielding approximation in their simula- 
tions, where a UV field is applied homogeneously to gas but 
with a cutoff at an assumed self-shielding density threshold 
of 10 -2 cm -3 . We have run an analogue to our HI simu- 
lation using the same self-shielding approximation instead 
of radiative transfer. Fig. [5] shows the neutral fraction ver- 
sus density phase diagram at z = 3 in this simulation. Apart 
from a much more discrete jump from ionized to neutral, the 
diagram is similar to the RHD counterpart (Fig. [4] bottom 
left), and we find 50% neutral fraction at half the density of 
the RHD counterpart, or at 0.025 cm -3 . In terms of getting 
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Figure 4. Neutral hydrogen at redshift 3 (HI, H2, H3 halos from left to right). Top row: Maps of projected Hi column density. Contours 
correspond to the lower limits for DLAs (N H = 2 10 20 cm -2 ), SLLSs (10 19 cm -2 ) and LLSs (1.6 10 17 cm -2 ). Bottom row: Phase 
diagrams of neutral hydrogen fraction, xm = nm/nn, versus density. The color scale represents mass- weighted probability per xm — tih 
bin over the plotted xm range. The quantization-like horizontal lines at the bottom of all diagrams are due to numerical precision of the 
cell variable xhii m RamsesRT, which is roughly seven digits. 




Figure 5. Phase diagram of neutral fraction versus density at 
redshift 3, in a version of the HI simulation where a self-shielding 
approximation is used (instead of RHD) of applying full-strength 
UV background at densities below 10 -2 cm -3 and zero strength 
above. 



but one should be careful not to move it higher than that to 
avoid over-predicting Lya emissivities due to photo-heating 
and photo-fluorescence. The approximation inaccurately de- 
scribes UV attenuation in more massive halos, where much 
of the gas is thermally ionized and thus UV transparent. 
This doesn't matter however, since an absence/presence of 
the UV background in gas which is already so ionized has a 
negligible effect on its Lya emissivity (which is dictated by 
collisional ionization equilibrium (CIE) heating/cooling. 



right the ionization state of gas at redshift 3 it thus ap- 
pears that this non-RT self-shielding approximation holds 
fairly well. One might perhaps consider moving the self- 
shielding threshold a factor of two towards higher density, 
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4 PREDICTED Lya LUMINOSITIES 

4.1 Computing the gas Lya emission 

In most astrophysical contexts, an electron in the excited 
level 2P of the hydrogen atom will practically instantly relax 
to the ground state (IS) via the emission of a Lya photon. 
There are two channels to produce such excited atoms, and 
hence to produce Lya radiation^] 

Collisional: A collision with a free electron excites the 
H-atom, which may release a Lya photon when it relaxes 
back to the ground state. The collisional emissivity is ap- 
proximated with 



Scoll — CLya(T) U e Tim ^Ly 



(4) 



where n e and nm are number densities of electrons and neu- 
tral hydrogen, respectively, and C\, ya {T) is the rate of col- 
lisionally induced lS-to-2P level transitions. An expression 
for this rate is given by G10, fitting results from Callaway 
(1987). It is always less than the hydrogen co llisional 



et al. 



excitation cooling rate, A colh used in the code (from|Maselli| 



et al. 2003), since cooling also takes into account excitations 
to atomic states other than 2P (the most likely of which 
is the non-Lya releasing 2S state). The ratio of C^ya/A^Ja 
goes from 71% at 10 4 K to 57% at 5 10 4 K. 

Recombinative: A free electron combines with a pro- 
ton at any level (^ 2), and may cascade down to the 2P 
level. The recombinative Lya emissivity of this process is 
given by 



0.68 a H i(T) rte nun £L y 



(5) 



where the 0.68- factor is the average number of Lya photons 



produced p er case B recombination (from Osterbrock & Fer- 
2006J and a^T) is the case B recombinations rate, 



land 



i.e. counting recombinations to all levels except directly to 
the ground one. We use the expression from |Hui &: Gnedin| 
(1997). 



Unless otherwise specified, the Lya emissivities calcu- 
lated in this paper are: 



Scoll + S r 



(6) 



Figure [6] shows the collisional and recombinative Lya 
emissivities of gas at typical stream density, assuming the 
gas is UV exposed (thick curves) and self-shielded (thin). 
Also shown in dotted black curves is the neutral fraction of 
the gas, approximately extracted using a simplified model 
of hydrogen-only and PIE/CIE equilibrium. 

The plot illustrates that it is crucial to be consistent 
in following the gas state of (T, xm,r) in the simulation 
code, since independently changing one of those factors with- 
out considering the effect on the others can have a dra- 
matic effect on s. If, for example, self-shielding is assumed 



2 To be exhaustive, there is a third one, which is absorption of 
photons with energies in the range 10.2 - 13.6 eV, which will 
excite the electron to any level ^ 2, which will in turn cascade 
down and sometimes produce a Lya: photon. This process is likely 
sub-dominant in the regime that we are investigating (iFurlanettol 
|et al.||2005| |Kollmeier et al.| [2010), and requires Lya radiative 
transfer, which we postpone to a future paper. 




Figure 6. Lya: emissivity of gas at number density nn 



3 10" 



which is close to the lower limit for accretion 



streams in our simulations. Thick curves show gas exposed to 
a UV field with r = 6.1 10 -13 s _1 , corresponding to redshift 
3. Thin curves show UV shielded gas (F = s _1 ). Blue dashed 
curves show recombinative Lya: emissivity, red solid curves show 
collisional Lya: emissivity and black dotted curves show the ap- 
proximate neutral hydrogen fraction in the gas, assuming equilib- 
rium between photoionization, collisional ionization and recom- 
binations. 



in post-processing and the neutral fraction changed accord- 
ingly without considering the change in temperature, the 
Scoii estimate can increase by almost an order of magnitude. 
This is unphysical - what really happens if gas suddenly be- 
comes UV shielded is that the temperature drops somewhat 
due to lack of photo-heating, the end-result being a slightly 
lowered value of e. 

The accuracy of the (T, xm, T)-state is secondary to 
consistency, because if the code handles things properly, e 
should simply reflect the work put into the gas by the UV 
background and gas advection. In the limit that the UV en- 
ergy input is negligible compared to gravitational heating, 
accurate modelling of the UV background isn't really cru- 
cial in the context of Lya emissivity, and applying e.g. a 
sensible shielding approximation like the one discussed in 
Sec. |3.3| should be OK. This breaks down when UV photo- 
fluorescence becomes non-negligible. 



4.2 Intrinsic luminosities 

Fig. [7] top row, shows maps of the rest-frame Lya surface 
brightness S of the three targeted halos, which is calculated 
by integrating the Lya emissivity (Eq. [6| along the LOS. 
We don't take absorption or scattering into account: These 
factors would certainly diminish the brightest spots associ- 
ated with CGM regions, but we don't expect them to affect 
the more diffuse streams much (see discussion in Sec. 4.3). 
The surface brightness is concentrated around CGM regions 
in all three halos, with S ~ 10 40 — 10 41 ergs -1 kpc -2 , and 
the brightness in streams is typically lower by one or two 
orders of magnitude. 

The bottom row in Fig.|7|shows radially cumulative Lya 
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Figure 7. Top row: Rest-frame Lycn surface brightness maps of the targeted halos (HI, H2 and H3 simulations, from left to right). 
Bottom row: Radially (3D) cumulative Lycn luminosities of the same halos. The thick black curves show the total luminosity. The 
dotted lines denote at what radius 50% of the total luminosity is accounted for. The blue curves show how the total luminosity is split 
between collisional and recombinative channels, given in Eqs. ^] and [5] respectively. The red curves show how the total is split between 
the contributions of the UV background and gravitational dissipation, discussed in Sec. |5.2| The thin black curves show the contribution 
of sub-CGM density gas (i.e. mostly streams) to the total luminosities. 



luminosities for the halos, i.e. fraction of the total luminosity 
within a given radius (black solid curves). Streams and the 
diffuse medium consistently contribute about 60% of the 
total luminosity, as indicated by the thin black curves. 

It is evident from both the maps and plots that there 
is a trend of more extended emission with increasing halo 
mass. In the HI halo, half of the total luminosity comes from 
the central 16% of the virial radius (dotted lines), while in 
H2 this radius is 20% and 33% in H3. Partly this is because 
the more massive halos consist of increasing quantities of or- 
biting galaxies so the surface brightness is just following the 
increased spread of CGM regions, as can be seen from red 
dots of surface emissivity in the maps and from correspond- 
ing steps in cumulative surface brightness in the plots. That 
is not the whole story though: The streams become more 
efficient Lya emitters with increasing halo mass. 

As seen from the blue curves in the luminosity plots, 
electron-hydrogen collisions dominate the total luminosity, 
and recombinations are borderline negligible, as should be 
expected outside ISM regions. This dominance increases 
with halo mass, with recombinations contributing 10% to 
the total in the HI halo and only about 5% in H2 and H3. 
The red curves in Fig. [7] will be discussed in Sec. [5] 



In Fig.[8]we plot total halo luminosities versus halo mass 
(which is denned, as in Table [I] as the total mass of all dark 
matter and baryons within the virial radius). From each of 
the three simulations we extract all halos from within the 
zoom-in volume and integrate Eq. [6] over their virial radii, 
excluding ISM gas. The halos roughly line up into a power 
law indicated by a red solid line, with exponent 1.25. There is 
a systematic tendency for halos in more massive simulations 
to be more luminous for a given mass, which is presumably 
an environment effect since the cosmological over- density of 
the zoom-in regions increases between the HI, H2 and H3 
simulations respectively. 

Large thick symbols mark the three main halos targeted 
in our simulations. For those we also plot luminosities ex- 
cluding consecutive phases of gas. Blue symbols show lumi- 
nosities when ignoring the ISM and the CGM, and green 
symbols show what happens if we also ignore the stream 
densities. The CGM accounts for about 40% of the total 
luminosity in all three targeted halos (see also Fig. [7| and 
the streams account for most of the rest, or 50 — 60%, with 
sub-stream densities accounting for 8%, 4% and 2% in the 
targeted halos of HI, H2 and H3 respectively. 

As can be read directly from Eqs. Q and (J5|, the Lya 
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Figure 8. Total halo luminosities versus halo mass for the three 
simulations. The halos are extracted from the zoom-in volumes, 
and we exclude sub-halos. The three targeted halos are indicated 
by large thick symbols. For those halos we also show luminosities 
excluding different phases of the gas. The red solid line indicates 
a power-law with exponent 1.25. 




Figure 9. Phase diagram showing Lya emissivity of gas versus 
density at z=3 in the HI halo. Each density bin is Lyo; luminosity 
weighted independently. The grey dashed lines show power laws 
that approximately fit the data. The color scale represents Lya 
luminosity weighted probability per e — tih bin over the plotted 
e range. 



emissivity of gas in principle scales with density squared, 
though temperature and ionization state have their influence 
as well. Fig.[9]shows a luminosity weighted phase diagram of 
Lya emissivity of gas versus density in the HI halo (the H2 
and H3 diagrams are similar). Over-plotted on the diagram 
in dashed grey lines are two power laws that the gas emis- 
sivity approximately follows, with a knee between 10 -2 and 
10 _1 cm -3 . The knee roughly corresponds to where the gas 
becomes self-shielding and the change in slope is caused by 
the corresponding transition in temperature and ionization 
state. Below the knee the gas emissivity is split in two ridges 
with slightly different slopes. The upper one has power in- 
dex « 2.5 and is dominated by collisional emission (Eq. 
whereas the lower one has power index ^2.2 and is domi- 
nated by recombinations (Eq. |5|. The emissivity above the 
knee is completely dominated by collisions. The stronger 
than 2 power law below the knee stems from the increasing 
abundance of neutral atoms with density and a tempera- 
ture that tends towards peak Lya emissivity, whereas the 
less than 2 power law above it results from the decreasing 
relative abundance of electrons with density. 



4.3 Observational properties 

We now consider mock observations of our simulated ha- 
los. To produce those we first convert the rest-frame surface 
brightness maps in Fig. [7] to observed surface brightness /, 
using 



Sfa 



4n(l + z) 



(7) 



where f a is a cosmological transmission factor that accounts 
for absorption and scattering of Lya photons on the LOS 
from the object to the observer. We adopt in this paper a 
value of f a = 0.66 based on the work of Faucher-Giguere 
et al. (2008). (f a is only applied to mock observations and 



not to the intrinsic emissivity and luminosity, Figures [7] and 
|8j). To the result of Eq. [7] we then apply a Gaussian point 
spread function (PSF) with a 0.6 arcsec full width at half 
maximum (FWHM) to mimic atmospheric and instrumen- 
tal distortion, and assume a camera pixel size of 0.2 arcsec 



(Fig. 10 ). This corresponds to very good seeing conditions in 
state-of-the-art instruments. We present maps made with a 
PSF about twice as broad in Appendix[C] These are directly 
comparable to the observations of Mil. 

Unlike FG10, we don't model the scattering of Lya pho- 
tons in this work. These authors show that Lya transfer 
dominates the spectral shape of extended Lya emission, but 
hint that it has little effect on the morphology and extent. 
Their Fig. 8 shows this to be the case for a halo correspond- 
ing in mass to our HI halo, if only the Lya emissivity of 
gas is considered - though their Fig. 9 also illustrates that 
strong point-like sources can produce extended Lya struc- 
tures via scattering. We will assume here that scattering 
has little overall effect on our predicted morphologies and 
extents, in the case that these structures are already well 
extended, though we do expect that it will likely produce 
subtle changes in observable LAB areas - indeed Fig. 8 in 
FG10 shows that the inclusion of scattering can make some 
observable Lya structures narrow down and others widen 
out. We will include and investigate the effect of Lya scat- 
tering on spectral shapes, luminosities, and morphologies of 
our objects in a future paper. 



4-3.1 Observed blobiness of cold accretion streams 

Mock observations of our three targeted halos are shown in 
the top row of Fig. |10| The middle contour in the maps 
is set at I -is — 10" 



erg s " cm " arcsec 2 , roughly cor- 



responding to current observation limits (e.g. Mil, Erb 
|et al.| |2011[ see Appendix [C] for a more accurate com- 
parison), and the inner and outer contours correspond to 
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Figure 10. Top row: Mock images showing predicted observed surface brightness in the targeted halos of HI, H2 and H3, from left 
to right. The contours mark 10 -17 , 10 -18 and 10 -19 erg s _1 cm -2 arcsec -2 . The images where computed using an optimistic PSF of 
FWHM 0.6 arcsec. Bottom row: Stellar density maps for the same halos, illustrating that the bright spots of Lya: emission are centered 
on galaxies. Over-plotted are contours marking 10 -18 erg s -1 cm -2 arcsec -2 in observed surface brightness (the middle contours from 
the upper maps). 



ten times brighter and ten times dimmer, i.e. 10 1 and 
10 -19 ergs -1 cm -2 arcsec -2 . 

Assuming I -is as our instrumental sensitivity limit, the 
HI halo (top left) is a Lya emitter that is centered on a 
galaxy, circularly symmetric in shape, about 20 kpc in di- 
ameter and doesn't trace streams. Thus the HI halo is not 
a LAB. The total observed luminosity, i.e. / integrated over 
the area within /-is, is L Q bs = 6 10 41 ergs -1 . 

The H2 halo observation (top middle) differs dramat- 
ically from that of HI. At I -is we do see a borderline gi- 
ant LAB, asymmetric and about 100 kpc in length, and we 
can see the end of an accretion stream poking out to the 
north-west. The observed luminosity integrated above I -is 
is L obs = 2 10 43 ergs -1 . The H3 halo (top right) has ob- 
servable Lya emission all over the place, is about 200 kpc 
in diameter and very asymmetric. Its observable luminosity 
is Lobs = 10 44 ergs -1 . 

Provided there is nothing special about these halos, we 
can conclude that in general the cooling emission from halos 
with masses greater than a few times 10 12 Mq can produce 
giant LABs (> 100 kpc) at redshift 3, assuming current in- 
strument sensitivity limits. Qualitatively this compares well 
with Yang et al.|(|2010j), who find that at redshift 2.3, LABs 



should occupy halos > 10 13 Mq. Qualitatively again, the 
maps presented in Sec. [C] which mimic the observational 
conditions of Mil, show that the morphologies of our simu- 
lated LABs are very similar to those observed. 

Interestingly, we note that the LABs produced by cold 
accretion streams are naturally extended in the direction 
of the main large-scale filaments that they are connected to. 



This is particularly visible for H2 and H3 in Fig. 10 (see also 



Fig. 14), and lends support to the observational findings of 
Erb et all (12011). 



Another matter are those mysterious LABs which do 
not seem to be centered on observed galactic counterparts 



e.g.lSteidel et al.|2000||Weijmans et al.|2010|[Prescott et al 



2011 ). We are not able to reproduce this phenomenon in our 



simulations. The bottom row of Fig. |10| shows stellar densi- 
ties in our targeted halos, with the I -is sensitivity contour 
over-plotted. Clearly all the peaks of Lya brightness would 
have continuum counterparts in observations, unless these 
counterparts would for some reason be hidden from view. 
Such LABs are rare among rare events, though, and our 
three simulations have little statistical chance of reproduc- 
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Figure 11. Plot of mock observed areas within contours of 
1.4 10 -18 ergs -1 cm -2 arcsec -2 versus halo mass for zoomed ha- 
los in the three simulations. Each halo is observed from three di- 
rections. The targeted halos are marked with thick symbols. The 
points are enveloped by two power-laws shown in the plot. 



ing such oddities. A larger sample of simulations would be 
required to investigate this issue further. 



4-3.2 Size distribution of simulated LABs 

We shall now statistically compare our results with a cata- 
logue of 202 observed LABs from the surveys described in 
Mil (courtesy of Yuichi Matsuda and team). The aim here 
is to derive a cumulative LAB area function from our results 
and see how it compares with real data. 

We follow Mil by assuming z — 3.1 in Eq. [7| and ap- 
plying a PSF with FWHM=1.4 arcsec. We calculate the ob- 
served LAB area A of each halo within the zoom regions of 
our simulations by integrating its total area above the sur- 
face brightness limit / = 1.4 10 -18 ergs -1 cm -2 arcsec -2 . 
We 'observe' each halo in three directions (x, y, and z). In 
Fig. |11[ we plot the LAB areas as a function of halo mass. 
The large thick symbols correspond to our targeted halos 
HI, H2, and H3. The observed LAB area is a reasonably 
well-behaved function of halo mass, with more massive ha- 
los producing larger LABs, and the points can be bracketed 
by a couple of power laws of indexes 0.87 and 1.19 (see Fig. 

n). 

We now make the assumption that extended Lya emis- 
sion is an inherent property of dark matter halos and that 
the observed LAB properties are direct functions of halo 
mass. This assumption is substantiated by our results (Figs. 
[8] and 11 ). We thus convolve the power laws of Fig. [II] with 
the halo mass function at redshift 3.1 (taken from Sheth 
\k Tormen 1999) in order to produce the cumulative area 
function envelope shown in Fig. |12[ There, the black dia- 
monds represent actual observations for comparison: they 
are a rough estimate of the area function at redshift 3.1 



Figure 12. The shaded region represents the boundaries of our 
predicted LAB area function, derived from the power laws in 
Fig. |11| The black symbols mark a rough area function derived 
from a sample of 202 observed LABs from the survey of Mil. The 
horizontally and vertically line-filled regions represent similarly 
predicted area functions, but with gas densities of nn ^ 0.3 cm -3 
and nu ^0.1 cm -3 excluded, respectively. 



based on the 202 LABs of the Mil survey, derived by bin- 
ning the LABs by area and dividing the count by the total 
survey volume of 1.57 10 6 Mpc 3 (the error bars are Pois- 
sonian). The comparison between our predicted area func- 
tion and the observationally derived one is very satisfactory, 
although we systematically over-predict the function by a 
factor of 2-3. 

There may be several causes to this over-prediction. 
First, our derivation of the observed LAB area function is 
too simplified. For example, we do not take into account the 
shape of the narrow band filter, or any 1/Vmax corrections. 
This introduces systematic errors that could well be of about 
a factor two. Second, our mock observations are also simpli- 
fied, and do not include noise, which could possibly affect the 
measured area in a systematic way. Third, perhaps we have 
overshot in our choice of f a = 0.66. As noted by G10, cosmic 
extinction may be stronger than average for sources that re- 
side in over-dense regions, as LABs tend to do. Fourth, our 
prediction is based on only a few objects and to a lesser de- 
gree the same applies to the observation-derived function. 
Fifth, the predicted LAB areas are sensitive to the applied 
PSF smoothing, which may not be entirely consistent in all 
the 202 observed LABs. And finally, we may lack physics in 
our simulations that would drive down the LAB areas. For 
example, Lya scattering, if applied, could induce a slight 
spread in the predicted rest-frame Lya surface brightness, 
which could in some cases bring down both the observed 
area and luminosity within sensitivity ordained brightness 
contours. Also, metal- line cooling may drive down the Lya 
emissivity of gas by cooling it below 10 4 K. Furthermore, as 



shown by van de Voort et al. (2011b 



(2011), and van de Voort &; Schaye 



Faucher-Giguere et al.| 
(2011|, feedback driven 
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Figure 13. Surface brightness profiles of our three targeted halos 
(solid) compared to the profile of LABd05 from |Prescott et al.| 
( |2011| > (dashed) and an average surface brightness profile of 11 
LABs from |Steidel et al.| ( |20TT| (dot-dashed). The thin solid lines 
represent different random orientations (50 for each halo) and the 
thick solid ones are averages of the thin ones. 



winds can destroy cold accretion streams in the vicinity of 
galaxies, hence terminating their Lya emissivities. 

Our predicted area function is not very sensitive to 
the density threshold of gas applied throughout this paper, 
where we have excluded ISM densities (tih > 1 cm -3 ) in 
our analysis. To illustrate this, Fig. [12] also shows, with line 
filled regions, predicted area function envelopes that have 
been derived from our simulations via a convolution with the 
Sheth-Thormen halo mass function, but including only more 
diffuse gas, nu < 0.3 cm -3 (i.e. sub CGM densities) and 
71h < 0.1 cm -3 , for the horizontal and vertical line- fillings 
respectively. The prediction using the sub-CGM densities is 
close to the original prediction, which can be expected since 
these densities account for ~ 60% of the total luminosities 
of all three targeted halos (see Figs. [7] and Even using 
71h < 0.1 cm -3 gas only (which is comparable to the more 
conservative prescriptions used in FG10) still produces gi- 
ant LABs hosted by massive halos and gives an area function 
that is compatible to the observational data. This confirms 
that the extent of our simulated LABs is largely driven by 
low density cold streams. 

We have also compared our results to observations via a 
LAB luminosity function (rather than the area function just 
discussed). However, since LAB emissivity typically peaks 
around compact sources, and since we neither model the 
emission nor absorption coming from the compact ISM re- 
gions, such a comparison is less robust than using the area 
function which should be more or less dictated by the state 
of more diffuse gas on much larger scales. The luminosity 
comparison, which is discussed in detail in Appendix [D] is 
actually surprisingly good, but it is problematic to draw any 
conclusions from it because of the lack of modelling of com- 
pact regions. 



4-3.3 Surface brightness profiles 

Fig. [13] shows observed surface brightness profiles of our 
three targeted halos. Profiles are taken for each halo in 50 
planes of random orientation, represented by thin coloured 
lines, and then these are averaged into the thick coloured 
lines. The profiles are transformed from rest-frame to ob- 
served surface brightness via Eq. [7] assuming redshift 3, 
but no smoothing or cosmic extinction is applied, and as 
before Lya scattering is neglected. Over-plotted are sur- 
face brightness profiles from observations: The black dashed 
curve represents a Sersic fit to the observed profile of the 



giant LABd05 at redshift 2.656 from |Prescott et al.| ( |2011| ), 
which we have scaled to z=3. Notably, LABd05 doesn't have 
a galactic counterpart at, or even close to the peak of Lya 
emission, though it has 17 small galaxies substantially offset 
from the peak (by > 20 kpc). The black dot-dashed curve 
is an exponential disk fit to an average of 11 LAB profiles 
observed at z ~ 2 — 3, reported in Steidel et al. (2011 ), with 
no scaling applied. 

The profiles of H2 and H3 are similar in shape and mag- 
nitude to the observed profiles. Interestingly, each of those 
compares favourably to different observations, with the H2 
profile being similar to LABd05 and the H3 profile similar to 



the 11 LABs from Steidel et al. (2011 ). The comparison indi- 
cates that these observations fit well within the model of cold 
accretion powered LABs, but due to the very limited statis- 
tics of our simulations (i.e. one halo per mass bin of three), 
and different redshifts of the observed LABs, it is problem- 
atic to make quantitative deductions, e.g. about masses of 
the host halos of observed LABs. Rather than representing 
different halo masses, the different profile shapes (and to 
some degree their magnitudes) may just as well reflect the 
different morphologies one may find in galactic groups and 
clusters. 



Prescott et al. (2011) compare the LABd05 surface 



brightness profile with simulated profiles from G10 and 
FG10, and find that the simulations appear to fit very badly 
with reality, with the G10 profile being both too peaky at 
the center and too shallow at large radius, and the FG10 
profile being too weak and steep. Our HI profile actually 
agrees with the simulated profile from FG10 (their model 7, 



see Fig. 9 in |Prescott et al. |2011| ), and it seems to us that 
FG10 in fact don't pose any mismatch with the LABd05 
observation: The fault lies in Prescott et al. (2011) assum- 



ing that the surface brightness profile scales linearly with 
halo mass, which is not at all the case judging from our 
simulations (and to be fair, these authors admit that their 
assumption is probably not accurate). 

We admittedly don't provide large statistics here, but 
we can conclude that the surface brightness profiles pro- 
duced by our simulations do not disagree with LAB obser- 
vations, and at the same time we can argue that neither do 
the simulations of FG10. 



4.3.4 Implications for future observations 

Having demonstrated reasonable agreement between our 
simulations and LAB observations, we now wish to highlight 
a prediction from our work which is particularly relevant in 
the context of direct searches for IGM emission at high red- 
shifts. The outermost contours in the upper row of Fig. [lO] 
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Figure 14. Same as Fig. |10| but showing the large-scale Lya: 
map of the H3 halo and its environment. Thick inner (thin outer) 
contours mark / = 10 — 18 (10 -19 ) erg s —1 cm -2 arcsec -2 . Large- 
scale streams connecting massive halos and extending over several 
Mpc would be visible at 10 — 19 ergs -1 cm -2 arcsec -2 . 



mark Lya brightness at 10 -19 ergs -1 cm -2 arcsec -2 . At this 
limit, accretion streams start to show up even in the least 
massive halo, and in the more massive halos we would de- 
tect them unambiguously. The deepest observations to date 
are not quite there yet, but almost, and this is an exciting 
perspective. Perhaps even more exciting is the map shown 
where the thin (resp. thick) contours again mark 



14 



on Fig. 

the limit at 10 -19 (resp. 10 -18 ) erg s -1 cm - ^ arcsec - ^. This 
zoomed out view of our H3 halo shows that deep Lya obser- 
vations around massive halos may even reveal the large-scale 
filamentary structure of the IGM on scales of a few Mpc ! 

Although such observations are extremely challenging 
(if not plain impossible) today, upcoming instruments, such 
as MUSE ( |Bacon et al.|2006) or K-CWI (Martin et al.|2010| ), 
should greatly increase our chances of observing directly this 
source term of galaxy formation in a very near future. 



5 WHAT DRIVES THE Lya EMISSION? 

Since Lya scattering and stellar feedback are not included 
in our simulations, the only possible power sources of Lya 
emission in our results are gravity and the UV background. 
We now look at how gravitational heating contributes to 
the Lya emission along streams and attempt to quantify its 
efficiency. We consider the contribution of UV fluorescence 
and show how it is sub-dominant for typical values of the 
UV background. We conclude this section by discussing to 
what extent locally enhanced UV fluxes could boost Lya 
emission from cold streams. 



5.1 Gravitational efficiency 

Gravitational heating is generally viewed as a progressive 
release of gravitational potential energy that heats the gas 
along the cold streams (|Dijkstra Loeb|2Q09|. As long as 



this heating is not too fast it can be balanced by radiative 
cooling, and as long as the gas remains metal-poor and at 
temperatures ~ 10 4 K, Lya emission is the dominant cooling 
mechanism, meaning that the thermal energy is efficiently 
converted into Lya photons. 

Gravitational heating in cold streams can be 



parametrized by the gravitational efficiency f g 



the 



fraction of the change in gravitational potential that 
dissipates into thermal energy during in-fall. The rest is 
converted into bulk kinetic energy, increasing the speed of 
the gas. A value of f gra v = 1 thus means perfect conversion 
of potential into heating, implying constant in-fall speed, 
and fgrav = means that there is no conversion into 
thermal energy and the stream should be in free-fall. 



|Dijkstra fc Loeb (2009) derive an analytic model in which 
fgrav > 20% is required if Lya blobs are to be driven by 
gravitational heating in cool flows. 

For each stream we can distinguish in our simulated 
halos in a given output, we can extract the stream speed 
profile v s tr(r) by following its core from end to end. Using 
this and a corresponding free-fall profile Vff(r) for a body 
starting at a position and speed identical to the outer end 
of the stream, we can estimate f gra v with 



fgrav (f) — 



v 2 ff (r) - v 2 str (r) 



J ff 



(r) 



(8) 



where Vinit is the speed at the outer starting position. 

We calculate an approximate free-fall profile for the 
stream by assuming static state and spherical symmetry, 
integrating the free-fall speed from the starting position to- 
wards the halo center using 



dv ff (r) 



1 GM(< r) 



v ff( r ) 



air, 



(9) 



where r is radius, G is the gravitational constant and M(< 
r) is the total halo mass within r. 

In practice we divide the halo mass into radial bins 
ri, where increasing i corresponds to decreasing radius, and 
solve Eq. [5] by recursively computing 

1 GM(< r <+ i) 



Vff(r i+1 ) = v f f(n) + 



viji) 



n), (10) 



where M(< r^+i) is the mass measured within r^+i in the 
simulation, and the initial condition is the stream speed at 
the outer end, Vff(r ) = v in i t . 

In Fig. [15] we show phase diagrams for the three tar- 
geted halos of gas speed versus radius (normalized to R V ir), 
where we exclude all but gas at stream densities, so that 
the streams can stand out more clearly. In the smallest halo 
(HI) the streams pop out nicely, smooth and undisturbed 
basically over the whole radius range, though they do di- 
lute a bit at the central 10% of R V i r - In H2 we can still see 
streams, but they are much more disrupted, and not distin- 
guishable within the central 20% of R V i r . In H3 only a few 
streams can be distinguished in the outer 40% of R V i r , and 
in the central ~ 20% they are completely destroyed. 

For those streams we can clearly distinguish in the di- 
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Figure 15. Phase diagrams showing speed of stream gas as a function of radius within targeted halos HI, H2 and H3 from left to right. 
The color scale represents mass weighted probability per speed-radius bin. Only stream densities are included here so that the streams 
stand out in the diagrams and are not drowned in the more diffuse gas. Over-plotted in yellow are free-falling speed profiles. Plots above 
each phase diagram show gravitational heating efficiency in the clearest streams, compared to the free-fall profiles. 



agrams, we have plotted in yellow the corresponding free- 
fall profiles, using Eq. [lO] which show approximately the 
speeds that the streams would follow were they in free-fall. 
Qualitatively it can be seen that the streams are close to 
free-fall, though usually they lag a little behind the free- 
fall profile, and conversely on some occasions we even see 
streams that seem to accelerate faster than free-fall (due to 
sub- halos and/or the inaccuracy of assuming static state and 
spherical symmetry in our free-fall calculation). 

We plot our estimates of f gra v using Eq.[8]directly above 
each phase diagram. For two of the three streams we have 
extracted in the HI halo we get a fairly consistent estimate 
of fgrav °° 0.1 from the halo outskirts towards the central 
~ 15% of R V ir, whereas for the third (and more diffuse) 
stream we get a value which is two to three times higher. In 
the H2 halo things are much messier, and for those fragments 
of streams that we can extract we find a large scatter in 
fgrav, going from negative values to about 0.3 (the initial 
large values are a numerical noise due to resolution in the 
phase-space). Finally, in the H3 halo, we can only extract 
two streams at the outer edges of the halo, one of them 
showing fgrav ~ 0.1 and the other accelerating faster than 
our free-fall approximation. 

It appears that gravitational heating is inefficient if seen 
only as a smooth and steady process along unperturbed 
streams as in HI. However, heating and subsequent release 
of Lya photons seems to be more efficient when it involves 
disrupted and wiggly streams. This also appears reasonable, 
since the gas at the core of a straight and unperturbed 
stream can flow virtually unopposed towards the central 
galaxy whereas if the streams are wiggly and disrupted there 
should be greater opposition from the surrounding hot and 
diffuse gas. 

It remains to be seen how much photo-fluorescence 
from the UV background is contributing the Lya emis- 
sivities compared to gravitational processes, both smooth 
and messy. Before comparing these factors we describe how 
they're derived from the simulation output. 



5.2 Computing the Lya contributions 

Since we store the photon flux in each cell we can easily 
keep track of the photo-heating and photoionization rates 
in the gas. If we assume that every photoionization leads to 
a recombination and that all the energy provided by photo- 
heating is released via collisional excitations^] we can esti- 
mate the UV contribution to the Lya emissivity in each cell 
as: 

£uv = 0.7 + 0.68 T ram e hya , (11) 

where 1~L 1 is the photo-heating rate, the 0.7-factor is the 
conversion efficiency of cooling into Lya photon^] and the 
second term on the right is akin to Eq. [5] We refer to Ap- 
pendix [B] for how to calculate the photo-heating rate. The 
UV contribution tends to be overestimated and can in fact 
be estimated higher than e in hot regions where collisional 
excitation is not the dominant cooling channel, but since 
these regions are Lya dim anyway this isn't a concern. 

The only other driver of Lya emission in our simula- 
tions are hydrodynamical processes which we can coin grav- 
itational heating. Thus the approximate gravitational con- 
tribution to Lya emissivity can be calculated in each gas 
cell as: 

Sgrav = max (0, e - £uv). (12) 

5.3 Gravitational heating vs. UV fluorescence 

Applying Eq. [TT] we calculate the UV contribution to Lya 
luminosity in each gas cell. We find that the relative UV con- 

3 The timescale for recombinations in streams is on the order of 
10 5 - 10 6 years, which is short compared to the timescale for the 
in-fall of streams in these halos, ~ 100 million years. The cooling 
timescale in streams is typically on the order of 10 4 to 10 5 years. 

4 This factor (roughly) represents the ratio of Cl vol {T) to the 
hydrogen collisional excitation cooling rate discussed in Sec. |4.1| 
It means that we assume 70% of the energy dissipated via cooling 
to go into Lya photons. 
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Figure 16. Top row: Fractional UV background contribution to the gas Lya emissivity in the three main halos, HI, H2 and H3, 
from left to right. Shown are mass weighted averages along the LOS, and everything below stream densities (0.02 cm -3 ) is ignored. 
Bottom row: Density distribution of the total Lya luminosity of the same halos, split into the UV (red) and the gravitational (blue) 
contributions. Green columns denote the Lyo: luminosity from maximum fluorescence. The shaded area represents stream densities. Note 
that the histograms have different scales on the y-axes. 




Figure 17. Ratio of the rest-frame Lyo; surface brightness assuming maximum fluorescence (i.e. infinite UV flux, »S max uv) to our fiducial 
model, for halos HI, H2 and H3, from left to right. Fluorescence considerably boosts the densest clumps at the price of dimming the 
most diffuse streams. 
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tribution becomes weaker with increasing halo mass, with 
the ratio to the total halo Lya luminosity going from 8% in 
HI to 2% in H2 and 1% in H3 (see bottom row of Fig. 
However, the relative UV contribution is generally stronger 
on the edges of the halos than near their centers. 

In the top row of Fig.[l6]we map the fractional UV con- 
tribution to the Lya emissivity in streams and CGM gas, 
and in the bottom row of the same figure we plot the den- 
sity distribution of the total luminosity, split into the UV 
(red) and gravitational (blue) contributions. As the maps 
and histograms show, the UV contribution is negligible ev- 
erywhere except for the smooth and diffuse streams with 
TtH < 0.05 cm -3 in the HI halo and on the outskirts of the 
H2 halo. A comparison with the mock observations in Fig. |10| 
reveals that these diffuse streams where the UV background 
contribution is non-negligible are nowhere close to being ob- 
servable and all the observable emission is completely domi- 
nated by the gravitational contribution. The UV background 
contribution to extended Lya emission can thus safely be ig- 
nored, at least until the observational sensitivity increases 
by two orders of magnitude or so. 

Inclusion of local stellar UV radiation in our simula- 
tions may boost the UV contribution, and thus both the 
total luminosity of the halos and the extent of observable 
emission. Alternatively, the presence of a luminous quasar 
nearby may also significantly enhance the Lya luminosity 
through fluorescence, as demonstrated by |Cantalupo et al.| 
(20051 andlKollmeier et al.l (120101). 



significantly to LABs, as demonstrated by |Cantalupo et al.| 
(2005} and Kollmeier et al. (2010[). However, from the per- 



5.4 Maximum fluorescence 

Even though we lack in this work the inclusion of local 
sources of UV radiation, we can still evaluate the upper limit 
to the fluorescent Lya luminosity that we can get from our 
simulated structures. This gives us a idea of the relative lu- 
minosity increase a local UV enhancement would provide, 
both in terms of the global luminosities of our halos and in 
terms of where the Lya emissivity is boosted and where it 
is dimmed, compared to the gravitationally driven emission 
we have calculated. 

To show the maximum fluorescent Lya luminosities we 
can obtain, we re-calculate the Lya emissivity of the simu- 
lated gas in the limit that xhii = 1 everywhere, correspond- 
ing to an infinite flux of UV photons. Note that in this limit 
Scoii is zero everywhere and the Lya emissivity is purely 
recombinative. 

The green columns in the histograms in the bottom row 
of Fig. [16] show the Lya luminosity of gas in different den- 
sity bins in this limit. In all three halos at nu > 0.1cm -3 , 
maximum fluorescence outshines the normal gas due to the 
increase in HII abundance, while at lower densities it is dim- 
mer than what we predict with gravitational heating, be- 
cause collisional emission goes to zero. 

In Fig. |17| we show the effect of maximum fluorescence 
on the Lya emissivity of the gas within our three simulated 
halos. This is displayed as the ratio of ^maxuv - the lu- 
minosity computed assuming an infinite amount of ionizing 
photons everywhere to S - the luminosity used in the 
rest of the paper, which assumes a standard (though inho- 
mogeneous) background value. Clearly, a strongly enhanced 
UV fluorescence will boost the Lya emission in a signifi- 
cant fraction of the gas (the blue part), and may contribute 



spective of observing accretion streams, the price to pay is 
the strong dimming of lower density structures (red) . 

This maximum fluorescence scenario is obviously opti- 
mistic, and only a tiny volume fraction of the Universe will 
likely come close to it, in the vicinity of rare and bright 
quasars in over-dense regions. Most of the IGM will more 
likely be in a regime comparable to our fiducial description, 
and its Lya luminosity will be powered by collisional exci- 
tation. 



6 SUMMARY AND CONCLUSIONS 

We have in this work addressed the questions of whether 
gravitational heating may be the main driver of LABs, and 
how close we are to making direct and unambiguous detec- 
tions of cold accretion streams via their Lya emission. 

To this purpose we have run and analyzed cosmological 
RHD simulations specifically tailored to accurately predict 
Lya emission from extended structures. These simulations 
are idealized in the sense that the effects of stellar feedback 
and Lya scattering are ignored such as to isolate the effi- 
ciency of gravitational heating in generating Lya photons. 
Our analysis is focused on redshift 3, which corresponds to 
most LAB observations. 

Our approach improves upon previous works in the fol- 
lowing ways: (a) Using RamsesRT, our newly developed RHD 
version of the Ramses code, we include on-the-fly propaga- 
tion of UV photons, which allows us to consistently and ac- 
curately model the self-shielding state in accretion streams 
and their resulting temperatures and ionization fractions, 
which are all very important to accurately predict their Lya 
emissivity. (b) We apply a novel refinement strategy that al- 
lows us to optimally resolve accretion streams to an unprece- 
dented degree and on much larger scales than previously. 
This allows us to spatially resolve the competition between 
gravitational heating and radiative cooling in those streams, 
(c) We post-process our simulation outputs with very small 
timesteps to ensure we also resolve said competition tempo- 
rally. Failing to do this leads to a dramatic underestimate of 
Lya emissivity of gas due to the commonly utilized numer- 
ical method of operator splitting, and previous works may 
have been marked by this problem, (d) We simulate more 
massive halos than hitherto done, based on the growing con- 
sensus that LABs are hosted by the most massive halos in 
the Universe. 

There are nontheless issues regarding uncertainties that 
potentially affect our results. One is the likely presence of 
artificial overcooling in shocks. As pointed out by |Creasey| 
et al. (2011), shocks - or the mean-free paths of particles 
inside them - are almost exclusively under-resolved in cos- 



mological simulations. The artificially broadened shocks can 
prevent the creation of hot and diffuse gas phase and in- 
stead allow it to efficiently cool and remain at temperatures 
where Lya emission is the most effective cooling channel. 
We may thus over-predict the Lya emissivities of shock re- 
gions in our simulations. However, this effect should be most 
severe in regions where gas is shocking on to galactic disks, 
and should thus be mainly constrained to CGM regions, and 
to the densest gas under consideration, i.e. nn ~ 1 cm -3 . 
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Weaker shocks may also exist at the boundaries of the dis- 
rupted streams in and around our more massive halos, but 
it seems unlikely that numerical overcooling is a big issue 
here, due to the high resolution, large volumes, low densi- 
ties, and the fact that the Lya emissivity is not particularly 
concentrated at the stream boundaries. 

It is an unavoidable fact that the denser the gas in our 
simulations, or in any simulations for that matter, the larger 
the uncertainty in its Lya emissivity. In particular, at den- 
sities > 0.1 cm -3 , gas may cool down to << 10 4 K via 
molecular or metal-line cooling, neither of which is included 
in our simulations. The densest gas is also in general the 
most Lya luminous in our simulations: What we term CGM 
gas (riH > 0.3 cm -3 ) consistently accounts for 40% of the 
total Lya luminosities of our halos, so we can estimate the 
total Lya luminosities to be uncertain by (very) roughly 
50%, and even more if we exclude still more diffuse gas than 
the CGM. We have however shown that our results and con- 
clusions regarding LAB areas are not sensitive to the density 
threshold applied (i.e. above which densities we ignore Lya 
emissivity) . 

Our main results are the following: 

• Cold accretion streams in halos more massive than 
~ 10 12 produces extended and luminous Lya nebu- 
lae which are by large compatible with LABs observed at 
z ^ 3, in terms of morphology, luminosity and extent. Grav- 
ity alone provides most of the energy, and we find that extra 
sources such as UV fluorescence, Lya scattering or super- 
winds are not necessary. This clearly doesn't rule out these 
other processes though, as they are likely all significant in 
the case of LABs, and further work is needed to study their 
complex interplay. 

• In our simulations, LAB area and luminosity are rea- 
sonably well-behaved functions of halo mass. We use these 
relations to compute the cumulative luminosity and area dis- 
tributions, and find that they are in reasonable agreement 
with observations given the relatively large uncertainties. 
This comparison however suggests that the combined effects 
of SN feedback, Lya scattering and an enhanced local UV 
field may possibly have a negative impact on the luminos- 
ity and extent of simulated LABs, when conjoined with cold 
accretion. 

• The model of gravitational heating as a driver of ex- 
tended Lya emission works according to our results, but we 
need to alter our notion of how it works: It is inefficient in 
the classic sense where gas accretion is smooth. Rather the 
accretion is messy and disrupted in massive halos and prob- 
ably involves some mass loss to the surrounding hot diffuse 
medium. 

• Our examination of maximum photofluorescence hints 
that in extreme cases local UV enhancement, e.g. near 
quasars, can boost the Lya luminosity of LABs and to a 
lesser degree their extent. As demonstrated by Cantalupo 
et al.| ( |2()05| and Kollmeier et al. (2010), this means that 
large accretion flows may be more easily observed in the 
proximity of quasars than elsewhere. 

• We find that cold accretion streams should be unam- 
biguously observable via direct Lya emission for the first 
time in the near future, on upcoming instruments such as 
MUSE and K-CWI which will allow to probe emission at 
surface brightnesses as low as ~ 10 _19 ergs _1 cm -2 arcsec -2 . 



Although we have significantly improved on previous 
work, a large number of theoretical issues remain to be ad- 
dressed. In forthcoming papers, we plan to investigate the 
effects of Lya scattering SNe-driven winds and local UV en- 
hancement from star formation. 
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Figure Al. Evolution with redshift of the Ly man-continuum 
part of the UV spectrum of |Faucher-Giguere et al.] ( |2009| which 
we use in our simulations. The plot shows photon flux versus 
wavelength (A) for selected redshifts. The vertical lines indicate 
how we split the spectrum into three (Hi-, Hel-, Hell-ionizing) 
photon packages. 



J(y) dv. 



(A3) 



When injected this way, the photons flow into adjacent 
cells which are above the UV density threshold and thus 
evolve into local photon fluxes representing the quasi- 
homogeneous UV field. Fig. |A2| shows how the package prop- 
erties evolve with redshift. 



APPENDIX B: CALCULATING THE 
PHOTOIONIZATION AND PHOTO-HEATING 
RATES 

The hydrogen photoionization rate F for the Lyman- 
continuum, in units of ionization events per hydrogen atom 
per unit time, is given by 



Jo 



crm(y)J(y) dv, 



(Bl) 



where am is the hydrogen ionization cross section and J is 
the local photon flux, integrated over all directions. Since 
the UV spectrum in our RHD simulations is discretized into 
three photon packages, the photoionization rate is extracted 
from each gas cell in the simulation output as 



APPENDIX A: THE QUASI-HOMOGENEOUS 
UV BACKGROUND 



We use the UV background model of Faucher-Gig uere et al.| 
(2009), which is available on the web, and consists of the 
redshift-evolving spectrum shown in Fig. |A1| As indicated 



by vertical lines in the plot, we discretize the spectrum into 
three photon packages; Hi ionizing with frequencies in the 
range (Vm, ^Hei); Hel ionizing in the range (vuei, ^Heii); and 
Hen ionizing in the range (z/Heii, oo). All photons belonging 
to a package i share the common properties of flux Ji , aver- 
age cross sections &ij , where j stands for the three ionizable 
primordial species Hi, Hel and Hell, and average energies e~ij 
per photoionization (again versus the three species). These 
properties are integrated from the redshift-dependent UV 
spectrum and updated every coarse timestep in the simula- 
tion. They are derived in the following way: 

For each package i that is denned for the fre- 
quency interval (vi ,vn) and gi ven the UV spectrum 
J(y) [photons cm -2 s _1 Hz -1 ] (Fig. Al ), we assign an aver- 
age photoionization cross section against each species j (Hi, 
Hel, Hen) as 



% p x J(y)dv ' 



(Al) 



where we use the expressions for o~j(v) from Hui & Gnedin 



(1997). Similarly, we assign to each photon package aver- 
age photon energies per photoionization event against each 
species: 



f" % * hv (jj(y)J(y) dv 
f"^ aj(v)J(v) dv ' 



(A2) 



where h is Planck's constant. The flux injected isotropically 
into each diffuse gas cell is derived for each package as 



r = E^ : 



(B2) 



where ami is the average hydrogen ionization cross section 
for package i and is the local flux of package i photons 
(see Appendix [A]). 

The photo- heating rate % 1 for the Lyman-continuum, 
in units of energy per time per volume, is given by 

HI, Hel, Hell 



E 



aj^F 1 (v) [hv — ej] dv, 



(B3) 



where we sum the photo-heating rates over the primordial 
ion species Hi, Hel and Hell. Here rij and is the number 
density of ion species j, o~j(v) is the species' cross-section, 
F 1 ' (y) is local photon flux, h is Planck's constant and €j the 
photoionization-threshold energies for species j. 

With the discretization of the UV spectrum into three 
photon packages (see Appendix [A| the integral in eq. (B3) 
becomes a sum: 



HI, Hel, Hell 



E 



3 

E 



(B4) 



where F? is the local flux of photons in package 2, and aji 
and eij are the photon package properties defined in ap- 
pendix [A] We plot the redshift evolution of the per-species 
ionization- and heating rates in Fig. A2 (right). 



APPENDIX C: MOCK LAB MAPS 

Fig. |Bl| shows mock observation thumbnails of the most lu- 
minous halos in the H2 and H3 simulations, produced in 
the same way as the ones in Fig. [To] but applying obser- 
vational parameters to match the surveys of Mil for di- 
rect comparison (these images can also be compared with 
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Figure A2. Left.Redshift-dependent properties of the three photon packages, integrated from the spectra in Fig. |A1| Right: The 
redshift evolution of the photoionization rates against each species (upper plot) and per-species photo-heating rates (lower). 



thumbnails in e.g. Yang et aL|2010 and Erb et aL||2011 ). In 
practice, this means that we assume our objects are at red- 
shift 3.1, smooth the images with a PSF with FWHM=1.4 
arcsec, and put the thick inner surface brightness contours 
at I = 1.4 10 -18 ergs -1 cm -2 arcsec -2 (the thin outer ones 
are at I = 10 -19 ergs -1 cm -2 arcsec -2 ). Morphologically 
our mocks resemble real LABs, asymmetric with a slight ten- 
dency to be filamentary and often having short sub-filaments 
that poke out of the main structure. In other words, our sim- 
ulated LABs look like real ones. 



APPENDIX D: COMPARING THE LAB 
LUMINOSITY FUNCTION WITH 
OBSERVATIONS 



In Sec. |4X2) we derive a cumulative area function from 
our mock LAB observations and compare to observations. 
We have done the same comparison for a cumulative lu- 
minosity function, though our prediction should be less ro- 
bust than the area function due to our lack of ISM mod- 
elling. Following Mil, we assume z = 3.1 and apply a 
PSF with FWHM=1.4. We calculate the observed Lya lu- 
minosities Lobs of all halos within the zoom regions of 
our simulations by integrating the surface brightness within 
/ = 1.4 10 -18 erg s -1 cm -2 arcsec -2 contours. These lumi- 



nosities are plotted against halo mass in Fig. Dl (left). Note 
the difference between L Q bs in this plot and L in Fig. [8] 
There we plot intrinsic luminosities of halos whereas here 
we plot observable luminosities, assuming instrument sen- 



sitivity and cosmological extinction (f a = 0.66). The plot 
points are enveloped by a shaded region bordered by two 
power laws, with indexes 1 and 1.45 as indicated in the plot. 
These power laws are convolved with a Sheth-Tormen halo 
mass function flSheth fe Tormen|1 999) at redshift 3.1 m or- 
der to produce the cumulative luminosity function envelope 



in Fig. Dl (right) 



The black diamonds in Fig. Dl (right) represent the 
observations from Mil, derived by binning the LABs by lu- 
minosity and dividing the count by the total survey volume, 
with Poissonian error bars. The comparison of our results 
to the observations is good, somewhat surprisingly so con- 
sidering the lack of modelling of the emission and absorp- 
tion in the compact peaks of Lya emission, which contribute 
substantially to the total luminosity. As with the area func- 
tion, we over-predict LAB abundances, though the predic- 
tion here is slightly closer to observations than in the case 
of areas. 

The plot also shows, with line filled regions, predicted 
luminosity function envelopes, where gas of densities nn > 
0.3 cm -3 and nu ^ 0.1 cm -3 is excluded from the analy- 
sis, for the horizontal and vertical line-fillings respectively. 



Much as with the area function (Fig. 12), excluding CGM 
densities and above (tih ^0.3 cm -3 ) has relatively small 
impact on the luminosity function. However using the lower 
density threshold of excluding gas with nu ^ 0.1 cm -3 gives 
a very abrupt change in the function, illustrating that the 
simulated LAB luminosities are more sensitive to the applied 
density cut than their areas. 
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Figure Bl. Mock observations showing / for the largest objects in the H2 and H3 simulations at redshift 3, smoothed with a Gaus- 
sian PSF of FWHM=1.4 arcsec to match recent observations. The physical widths of these thumbnail squares are all identical at 
300 kpc (« 40 arcsec). The thick inner contours correspond to I = 1.4 10 -18 erg s -1 cm -2 arcsec -2 and the thin outer ones to 
10 — 19 ergs -1 cm -2 arcsec -2 . 



APPENDIX E: COMPARISON TO PREVIOUS 
THEORETICAL WORK 

Our work is similar in nature to the work of FG10 and G10 
(see Sec.[T]), and for comparison to those we have chosen the 
mass of the HI halo to be similar to the halos on which they 
focus their analysis. 

Fig. 12 in G10 shows a plot of halo luminosities ver- 
sus virial masses at redshift 3.1, much like our Fig. [8] Their 
mass- luminosity power-law exponent is ~ 0.8, which is con- 
siderably shallower than our value of ~ 1.25. Their less 
massive halos are more luminous than ours, with our ones 
catching up around 10 12 Mq. Their targeted halos of masses 
« 4 10 n M o are typically a few times more luminous than 
ours, and much more extended in Lya emission. Their Fig. 
10 shows mock observations of two of their targeted halos. 
Also using /-is as a sensitivity limit, they have observable 
Lya emission which is very asymmetric, clearly traces ac- 
cretion streams and extends to about 100 kpc in length. 

Their prediction is probably a bit over the top, since a 
giant LAB in a halo of this size implies that LABs should 
be very common, and this contradicts the generally accepted 



view that they are uncommon and associated with unusually 
over-dense regions in the Universe (e.g . |Steidel et a"T1|2000| 



Prescott et~aT1|2008| |Yang et al.1|2010| ). The cause of their 



over-prediction can probably be traced to an overestimate in 
gas temperatures due to their self-shielding approximation 
being applied in post-processing, as pointed out by FG10 
(see also our Sec. 



4.1). 



Fig. 2 in FG10 shows a plot of halo luminosities ver- 
sus halo masses at redshift 3 for the various numerical ap- 
proaches. Their estimate which is most comparable to ours 
(their prescription 7, that sums all gas) has a power law 
with exponent ~ 1.1, steeper than that of G10, but still a bit 
shallower than our exponent, and their halos are slightly Lya 
dimmer than ours, typically around half the luminosity for a 
given mass, though this varies quite a lot due to scatter. The 
luminosity difference may be partly explained by the over- 
density of our simulated regions, which tends to increase the 
brighness of halos of similar mass, going from the least mas- 
sive to the most massive simulation. Their more conservative 
prescription excludes star-forming gas from their analysis, 
which in their simulations is gas with 71h > 0.13 cm -3 , likely 
a more realistic threshold than our nu > 1 cm -3 . On this 
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M vlr [l\y L obs [erg s- 1 ] 

Figure Dl. Left: Mock observed Lya luminosities within contours of 1.4 10 -18 ergs -1 cm -2 arcsec -2 versus halo mass for all zoomed 
halos in the three simulations. Each halo is observed along three simulation box axes and thus is represented by three points in the plot. 
The targeted halos are marked with thick symbols. The points are enveloped by two power laws shown in the plot. Right: Predicted 
luminosity function of LABs (shaded region), according to the power-law envelope from the left plot. The black symbols mark a rough 
luminosity function derived from a sample of 202 observed LABs from the survey of Mil. The horizontally and vertically line- filled 
regions represent predicted luminosity functions with gas densities of nn ^ 0.3 cm -3 and tih ^0.1 cm -3 excluded, respectively. 



exclusion, the luminosity drops by 1-2 orders of magnitude. 
This is a bit more dimming than we find in our results: If the 
left plot of Fig. |Dl| is considered, where the region filled with 
vertical grey lines corresponds to our luminosities where gas 
with 71h > 0.1 cm -3 is excluded from our analysis, it can 
be seen that the luminosity drops by < 1 order of magni- 
tude compared to our normal prescription of including all 
gas with riH > 1 cm -3 . 

They also show mock observations of a 2.5 10 11 Mq 
halo at redshift 3, that includes Lya scattering. Their Fig. 
7, middle left, can be compared to ours (again, their pre- 
scription 7, that sums all gas). A contour at I -is marks a 
very circular source centered on a galaxy, about 15 kpc in 
diameter, which is similar to our HI halo observation. The 
Lya luminosity of their halo is 8 10 41 ergs -1 , close to the 
'observed' luminosity of our HI halo of 6 10 41 ergs -1 . 

In terms of the emission coming from ~ 10 11 Mq halos 
we thus seem to be in fair agreement with FG10, though the 
LABs produced by our simulations appear to be somewhat 
more luminous, even when matching their more conservative 
prescriptions. In terms of LAB extent it is harder to tell, 
since they don't consider mock observations of halos more 
massive than 3 10 11 Mq, and it is thus hard to tell whether 
or not their massive halos produce LAB-like objects. 
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